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ABSTRACT 

We present results from a pulse phase resolved spectroscopy of the complex 
emission lines around 1 keV in the unique accretion powered X-ray pulsar 
4U 1626-67 using observation made with the XMM-Newton in 2003. In this 
source, the red and blue shifted emission lines and the line widths measured 
earlier with Chandra suggest their accretion disk origin. Another possible 
signature of lines produced in accretion disk can be a modulation of the line 
strength with pulse phase. We found the line fluxes to have pulse phase depen¬ 
dence, making 4U 1626-67 only the second pulsar after Her X-l to show such 
variability. The O VII line at 0.568 keV from 4U 1626-67 varied by a factor of 
~4, stronger than the continuum variability, that support their accretion disk 
origin. The line flux variability may appear due to variable illumination of the 
accretion disk by the pulsar or more likely, a warp like structure in the accre¬ 
tion disk. We also discuss some further possible diagnostics of the accretion 
disk in 4U 1626-67 with pulse phase resolved emission line spectroscopy. 

Key words: X-ray: Neutron Stars - accretion, accretion disk - pulsars, indi¬ 
vidual: 4U 1626-67 


1 INTRODUCTION 


4U 1626-67, discovered with the Uhuru satelli te (Giacconi e 


pulsar with a pulse period of 7.7 seconds flRappaport et al. 


al 


197 


1972) is an accreting X-ray 
7). It is one of the very few 
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X-ray pulsars in low mass X-ray binary (LMXB) systems. The magnetic fie 


neutron star is estimated to be of B rj 4x 10 12 Gauss (j Coburn et al. 


2002 


d str ength of the 


Orlandini et al. 


1998 


Iwakiri et al 


20121 ) from the cyclotron line features in its X-ray spectrum at around 
37 keV. This system with small mass function of / < 1.3 x 1O~ 6 M0 h as a neutron star 


that accretes from an extremely low mass compan i on (0 0.1 Mq) (IMiddleditch et al 


1981 


Levine et al. 


1988 


Chakrabartv 


1998 ). 


Jain et al. 


(120081 ) reported an upper limit of 10 lt- 


rns for pulse arrival delay using RXTE data. The optical light curve shows reprocessed 


X-ray pulses at 130.38 mHz with multiple sideba nds. The sic 


revealed an orbita’ 


period of 42 minute ( IMiddleditch et al. 


ebands in the power spectrum 


1981) that was confirmed by 


Chakrabartv! (119981) . The shor t orbital period suggests that the binary system contains a 


hydrogen depleted companion dPaczvnski fc Sienkiewicall9811). The optica 


strong UV excess and high optical pulse fraction (McClintock et al 


1977|, 


counterpart has 


19801 ). 4U 1626- 


67 under went two torque reversals since it s discovery. It was initially observed in a spin- 


up state ((Levine et al. 


1988 


star began to sp in down (IWilson et al 


Nagase 


989f). this trend reversed i n 1990 and the neutron 


1993 


Bildsten et al. 


Chakrabart v et al. 


1 994 ll 9971 : 


Prince et al. 


1994; 


19971 ). After the steady spin-down phase of about 18 years, the neu¬ 


tron _star_again underwen 


(Ijain. Paul fe Dutta 


; a torque reversal anc 


2010 ; 


Camero-Arranz et al 


made a transition to spin-up state in 2008 


2 0101 ). The pulse profiles of this source show 


a strong correlation with its torque reversals. The pulse profiles during two s pin-up eras were 


quite similar but they were different in the spin-down era (lBeri et al. 


20141 ) 


The spectrum of 4U 1626-67 has been studied many times using data from c 


observatories. Observa 


ASCA (Align 


(Krauss et al. 


2012 ; 


ini et al. 


Iwakiri et al. 


ions with HE AO 1 (IPravdo et al. 


1995 ). Beppo-SAX (Orlandini et ahlll998l ). C handra , XMM-Newton 


2007), RX TE (iJain. Paul fc Dutta! l2Qlol and Suzaku (ICamero-Arranz et al. 


1979 ). Tenma (Kii et al 


Terent 


1986 ). 


2012 ) showed two component spectrum consisting of a blackbody and a 


power-law. These continuum spectral parameters showed variation during the eras of torque 
reversals. The blackbody temperate 
sition from spin-up to sp in-down p 


Aneelini et al. 

1995; 

Vaughan & Kitamoto 

1997 


Owens. Oosterbroek & Par mar 

1997 

), followed by transition to value of ~ (0.5-0.6) keV dur- 

ing second phase of spin-up ( 

Jain. Paul & Dutta 2010; Camero-Arran 

lz et al. 

2012 

). During 

the first spin-up phase the value of powerlaw photon index was 1.5 ( 

Vaughan & Kitamoto 


1997 ). The spectrum became relatively harder during the spin-down phase with the powerlaw 
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19971; 


ton index of 0.8-1.0 was measured fl.Tain. Paul fc Putt 


n index of 0.4-0.6 

(Angclini et al. 

1995 

; Vaughan & Kitamoto 

1997; 

Owens. Oosterbroek & Parmar 

Orlandini et al. 

1998; 

Yi & Vishniac 

1999 

) and after the second torque reversal, a pho- 


2010 


Carnero-Arranz et al. 


20121 ). 


The presence of unusually bright Neon (Ne 

and Oxvgen (O) lines were reportec 

from 

many spectroscopic observations 

(Angclini et al. 

1995 

Owens. Oosterbroek & Parmar 

1997; 

Schulz et al. 

2001; 

Krauss et al. 

2007 

). The equivalent 

width and intensity of these emission 

lines were found to be variable (( 

Jamero-Arranz et al. 

2012 

). Observations made with Chan- 


dra revealed a double peake d nat ure of these line features, indicating their formation in the 
accretion disc flSchulz et al.l 200li ). The presence of these emission lines is quite interesting 
as these are not observed in any other X-ray pulsar. Oxygen emission lines are observed 


in both UV and optical spectra. However, optica 


emission features flHomer et al. 


2002 


Werner et al. 


spect rum of this source lacks in Neon 


20061 1. 


There was a gra dual 


2005 (Krauss et al 


decrease in the X-ray flux in the 0.3-10.0 keV band during 1977- 


2007), followed by an incr ease in 


sal (IC ame ro-Arr an z et al. 


2010 ; 


Jain, Paul & Dutta 


he flux during the second torque rever- 


2010! ). Observations made during spin- 
down period indicated a variations in flux of low energy emission li nes w ithin the same phase 


of torque reversal. The line intensities measured using ASCA data (lAngelini et al.lll995l) wer e 


higher co mpared to that mea s ured using BepyoSA X flOwens. Oosterbroek fc Parmar 


Chandra (Schulz et al 


2001 


Krauss et al. 


Suzaku ( Carnero-Arranz et al 


2007), XMM-Newton ( Krauss et ah 


997!), 


2007) and 


2012). The fluxes of these emission lines increased during 


transition from spin-down to spin-up (ICamero- Arranz et al .1120121 ). 


Pulse-Phase resolved spectroscopy is an important tool to probe into the geometry of the 


(Pravdo et al. 

1979; 

Kii et al. 

1986) 


19861) . A phase dependence of 3 7 keV Cyclotron Re sonance 


Scattering Feature (CRSF) was reported using Suzaku data by 


Iwakiri et ah 


( 20121 ). where 


they also observed it to be as an emission feature at certain phases. 


We present here a pulse phase resolved spectral analysis of 4U 1626-67 using XMM- 
Newton EPIC-pn observation made in 2003 to show the behaviour of the O VII (0.568 keV), 
O VIII (0.653 keV), Ne IX (0. 915 k eV), Ne X (1.021 keV) lines in different spin phases. Using 
ASCA data, lAngelin i et ah. (19951 ) reported that Ne Lya has no spin-phase dependence. 
Here we present more detailed analysis of the emission line fluxes in narrow phase bins using 
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XMM-Newton observation which has higher sensitivity than ASCA. The paper is organised 
as follows: The details of the observation and data reduction procedures are mentioned in 
section 2, section-3 includes timing analysis of the source, followed by the details of the 
spectroscopy. Finally, the implications of our results are discussed in section 5. 


2 XMM-NEWTON OBSERVATIONS AND DATA REDUCTION 


XMM-Newton operates in the energy range of 0.1-15 keV. The instruments on board XMM- 
Newton are European Photon Imaging Camera (EPIC), Refle ction Grating Spectrometer 
(RG S) and Optical Mon itor (OM). EPIC consists of two MOS ( Turner et al. 2001 1 and one 


PN (IStriider et al. 


200111 CCD arrays having moderate spectral resolution and a time res¬ 


olution in the range of micro-seconds to 2.6 seconds depending on the instrument and the 


mode of observation. High-resolution spectral i nformation 


Reflection Grating Spectrometer (RGS) (Iden Herder et ah 


E/dE ps 300) is provided by the 
20011 ). Simultaneous optical/UV 


observations are made with Optical Monitor (OM). 


We have used the longest observation (84 ks) of 4U 1626-67 made with XMM-Newton on 
20 August 2003 (ObsID-0152620101). The EPIC-pn data of this observation were collected 
in the small window mode and using the medium filter. 


HEASOFT 6.12 and SAS-12.0.1 were used for the reduction and extraction. For the 
extraction of the cleaned EPIC-pn events, standard hlters were applied. The observation 
was checked for the particle background using the SAS tool evselect and background flares 
were filtered out prior to spectral extraction. The event hies were then checked for pile-up, 
using SAS tool epatplot and as EPIC-pn has a fast readout time of 5.7 milliseconds for small 
window mode it was found that there is no problem of pile up in the source spectrum. 

For the extraction of the source spectra and lightcurves, circular region of 40 arcsec- 
ond radius was selected around the source centroid. For the background spectra and light 
curve extraction, two off-source regions were used with a radius of 40 arcseconds. Spectral 
extraction was done with PATTERN ^ 4 and FLAG=0. For reprocessing of the data and 
generation of the response hies updated version of the current calibration hies (ccf) [] was 
used. 


1 http://xmm2.esac.esa.int/external/xmm_sw_cal/calib/ccf.shtml 
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3 PULSE PROFILES 


The energy depenc 


Krauss et al. 


ent pulse profiles from this observation have already been reported in 


( 20071 ). As a prerequisite to the phase resolved spectrsocopy reported in the 


next section, we have repeated the period analysis. We have derived a pulse period of 7.67548 


seconds which is consistent with 


Krauss et al. 


(120071 ). However, the energy resolved pulse 


profiles shown in the left panel of Figure |T] in different energy bands show substantially more 
fine structures. 


Pulse profile of 4U 1626-67 is known to have strong energv dependence 

(Rappapor 

et al 


1977; 

Joss. Avni & Rappaport 1978; Pravdo et al. 1979 

Eisner, Ghosh & Lamb 1980; 

Kii et al. 

1986; 

Levine et al. 

1988 


Mihara 1995; i 

Angelini et al. 

1995; 

Krauss et al. 

2007 Jain et al 


2008; 

Iwakiri et al. 

2012 

)• 

Angelini et al. 

(1995); Krauss et al 

(2007) and 

lain et al. 

(2008) 


Jain. Paul & Duttc 

(2010 

) also reported pulse profiles similar to Figure Q] during the spin 



down phase of the pulsar 


two spin-up eras (IBeri et al 


h ese profiles were significantly different from that obtained in 


.2014). 


4 SPECTROSCOPY 

4.1 Pulse Phase averaged Spectroscopy 

Pulse-phase averaged spectral analysis was performed using the mean spectrum extracted 
from the EPIC-pn data, which was then grouped using ftool grppha with minimum 500 


counts/bin. The spectral fitting was done using Xspec Version-. 12.8.0 (lArnaudl 11996). The 
energy band chosen for fit was 0.3-10.0 keV. 


The continuum o 


and a power law (|Pravdo et al 


the phase averaged spectrum is well described by a soft thermal component 


1979 


Kii et al 


1986; 

Angelin 

et a 

1995; 

Owens. Ooste 

rbroek & Parmar 

; Krauss et al. 

2007; 

Jain. Paul & Dutta 

2010 



1997; 


Orlandini et al. 


Cainero-Arranz et ah 


1998: 


20121 : 


Schulz et al. 


20011 : 


Iwakiri et al 


2012 ). The two components cross over at 0.8 keV. 


We used bbodyrad model for fitting the soft part of the spectrum. 

Using only the continuum model shows significant excess in the form of emission line 
features around 0.5 keV and 1.0 keV (see Figure [2]). The RGS data from the same XMM- 


Newton observation was used by 


Krauss et al. 


(20 07 ). where they found four emission lines 


at 0.568, 0.653, 0.915, 1.02 keV. All the four lines were found to be narrow (a ~ leV). 
Following those results we added four Gaussian lines with energies fixed at 0.568 (O VII 
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Table 1. Best-fitting parameters of Phase averaged spectrum of 4U1626-67. 


Parameter 

Model Values 

N H (10 22 atoms cm -2 ) 

0.102 ±0.001 

Powlndex (r) 

0.793 ±0.003 

N a 

1 PL 

0.00656 ± .00004 

BBody (kT) keV 

0.229 ±0.001 

N b 

1 BB 

441 ± 18 

LineFlux 0 


O VII (0.568 keV) 

6.4 ±0.1 

O VIII (0.653 keV) 

3.7 ±0.1 

Ne IX (0.915 keV) 

1.36 ± 0.09 

Ne X (1.02 keV) 

1.88 ± 0.07 


Notes: Errors quoted are for the 68 % confidence range. 

a —y Power law normalisation ( Npj ,) is in units of photons cm -2 s —1 keV -1 at 1 keV 
b —> Blackbody normalisation (Nbb) is in units of photons cm -2 s —1 keV -1 
c —> Gaussian normalisation is in units of 10“ 4 photons cm -2 s —1 


He a ), 0.653 (O VIII Ly a ), 0.915 (Ne IX He a ) and 1.02 keV (Ne X Ly a ) to the continuum 
model, with their line widths fixed to <7=1 eV. On adding these lines an improved good 
fit was obtained with reduced y 2 (x 2 ) °f 1-175 for 1224 dof. The best fit parameters are 
given in Tabl e-1. These fit parameters are in good agreement with the values reported by 


Krauss et ah 


( 20071) using RGS data from the same observation. 


4.2 Phase-Resolved Spectroscopy 

The energy dependence of the pulse profiles in Figured! suggests a clear dependence of the 


energy spectrum on the sp in phase. A st r ong p 


67 has been r eport ed by 


Pravdo et ah 


Iwakiri et ah 


(119790 . 


r ase depend e nce o f t he continuum o 


Kii et ah 


(4986). 


Coburn et al. 


4U 1626- 


(120000 and 


(120120 using the phase resolved spectra acquired by HEAO-1, Tenma, RXTE 


and Suzaku respectively. 

In order to investigate the correlation between the features observed in the pulse pro¬ 
files with Neon and Oxygen emission lines in the spectrum of 4U 1626-67 we performed a 
pulse-phase resolved spectroscopy. For the phase resolved spectral analysis, the spectra were 
extracted with a phase bin of 0.05, using the ’phase’ filter of /tools fj task xselect. We 
used the same background spectrum and response hies as were used for the phase-averaged 
spectrum. X-ray spectra in different phases were grouped using ftools task grppha with 
minimum 250 counts per bin. For fitting the phase-resolved spectra we followed exactly the 
same technique as we opted for the phase-averaged spectrum. 


1 http://heasarc.gsfc.nasa.gov/ftools/ 
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The continuum parameters (Figure ]3j) show variation with the pulse phase. The equiv¬ 
alent hydrogen coulmn density (Nh) showed significantly higher value compared to the 
average value in the phase range 0.7-1.15 with some fluctuation. Temperature of blaekbody 
varied within the range of ~ 0.16-0.30 keV showing lower value around phase 0.7-1.15. The 
normalisation of blaekbody also exhibits a behaviour that is anticorrelated with its temper¬ 
ature profile. However, we note that these spectral parameters can have degeneracy between 
them and the rapid fluctuations of these parameter as seen in Figure [3] does indicate the 
presence of significant degeneracy. We therefore do not comment further on this. The pow- 
erlaw index however is determined from a broader energy range and shows lesser bin to bin 
fluctuation. We find a pulse phase dependence of the power-law index with spectrum being 
hardest near the pulse peak and it being soft near the dip of the pulse profile. 

Right panel in Figured] shows the behaviour of total flux, blaekbody flux, power law flux, 
Ne and O line intensities as function of the pulse phase. The total and power law flux show 
a sharp dip at around phase 1.0 along with two shallow dips near phase 0.3 and 0.7. This 
trend is quite similar to that of its pulse profile. Blaekbody flux also shows some structures 
in its profile with a dip between 0.2-0.3 phase. The behaviour of line fluxes shown in Figured! 
is compared to the features in its pulse profile. 

l.)The O Vll (0.568 keV) emission line shows strong variability by more than a factor 
of 4. This line is in maximum between the phase 0.0-0.15 and then decreases reaching a 
minimum intensity at around phase of 0.3. In the rest of the pulse phase also it shows some 
variation. It is interesting to note that peak intensity of this line is coincident with dip in 
the pulse profile. 

2. ) The O VIII (0.653 keV) shows some variation with spin phase but without any clear 
structures in it. 

3. ) Ne IX (0.915 keV) shows some variation with spin phase with a possible correlation with 
the pulse profile in 0.8-1.2 keV energy band. 

4. ) The Ne X (1.021 keV) line is consistent with no pulse phase dependence. 

The equivalent width (EW) of these emission lines are shown in Figure [U The EW of 
0.568 keV line also shows a maximum around phase 0.0-0.15 with a subsequent decreasing 
trend till phase 0.3. Some variation in the EW of this line is seen throughout the pulse 
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phase. The EW of 0.615 keV line also show variation with phase. A clip in the EW profile of 
0.915 keV emission line is coincident with the main dip observed in the pulse profile as well 
as its flux profile. However, the EW of 1.02 keV emission line is almost constant throughout 
its pulse phase. 

To evaluate the statistical significance of variation of the O VII (0.568 KeV) line with 
spin phase, we fixed the norm of the line in the phase resolved spectra with the value 
obtained from the phase-averaged spectrum. Fixing the line intensity to its phase averaged 
value resulted into an increase in the chi-squared (y 2 ) by upto 36 for 230 degrees of freedom 
(dof). In 5 of the 20 phase bins, the chance probability for the increase in y 2 is determined 
by F-test to be less than 10 -4 and in one of the bins it is as low as 5 x 10~ 8 , which clearly 
demonstrates the variation of the O VII emission line. 

A constant fitted to the flux (Figure CD) and equivalent width (Figure 0D of the O VII line 
returned y 2 values of 46 and 32 respectively for the 20 phase bins, which is another measure 
of the significance of the flux variability. 

To further inspect the significance of the line intensity variation, we have plotted the fraction 
of the number of phase bins for which the line flux differs from the average value by more 
than No - as a function of N. The same is compared with what is expected without any 
intrinsic variability of the line (Figure EJ. For any phase bin there is only 0.27 % chance 
probability for a 3 a deviation while in three of the 20 phase bins in Figure [5] the line flux 
deviates from the pulse averaged value by more than 3 a showing that the variability is 
statistically significant. 

We have also made a direct comparison of the spectra extracted in the pulse phase ranges 
0.0-0.15 and 0.15-0.4. At these phases, flux of this line shows maximum variation. In Figure[6l 
we have shown the best fitted spectra along with the model components in top panels for 
these two phase ranges. In the middle panels of the respective figures, we have shown the 
residuals obtained for the best fit and in the bottom panels the residuals are shown after 
setting the model line fluxes to zero. The residuals in the bottom panels clearly show the 
flux difference for the 0.568 keV line. 

To further investigate significance of the detected variation of the O VII line flux we 
simulated 10,000 spectra using the model parameters of the phase averaged spectrum, with 
an exposure time typical of a phase resolved spectra (2921 seconds). In each simulation, all 
the spectral parameters were taken to be same as their phase averaged values. Thereafter, 
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these simulated spectra were fitted using the same model allowing all but the energy and 
widths of the four lines to vary. A constant fitted to the flux of O VII line returned a y 2 
value of 8586 for 10,000 data points. We found that standard deviation of the O VII line 
flux variation in the simulated spectra was merely 15.2%, and 90% of the simulation results 
had flux within -36% and +40% of the average value. Thus the variation of line flux found in 
the simulated spectra is quite low compared to upto a factor of 4 variation observed in the 
O VII line flux, indicating that degeneracy between the spectral parameters did not result 
into values of the narrow line fluxes. 

Moreover, we have also performed similar analysis with different trial pulse periods of 
6, 7 and 8 seconds. The phase resolved spectra created using these trial periods were fitted 
using the same model. We did not find any significant variation of O VII line flux at these 
arbitrary trial periods. The constant fitted to the flux of O VII line flux shown in Figure |T| 
returned y 2 values of 16, 16 and 20 for 19 dof with the periods 6, 7 and 8 seconds respec¬ 
tively, compared to 46 for XMM-Newton data. This probe further strengthens our result of 
variation of O VII line flux accross the pulse phase. 

We also created pulse profiles in very narrow energy bands (of range of about 0.1 keV) 
to compare the shape of the profiles near 0.568 keV with that of the higher and lower energy 
continuum. However, because of the fact that the equivalent width (EW) of these emission 
lines is quite low (10-30 eV) in comparison to the energy resolution of about 100 eV and 
these line energies contribute to a small fraction of emission even for a band width of 100 eV, 
the difference between the pulse prohles was not very prominent, though we could observe 
some differences in the sharpness of the dips in the prohles. This was also consistent with 
our results obtained from pulse phase resolved spectroscopy. 


5 DISCUSSIONS AND CONCLUSIONS 


In the present work, we have performed a detailed pulse-phase resolved spectroscopy of the 
unique high magnetic held LMXB pulsar 4U 1626-67 using data from EPIC-pn onboard 
XMM-Newton observatory to investigate t he hux variation of complex emission lines arising 


from the highly ionised s pecies of Ne and O (lAngelini et al 


19971: 


Krauss et al. 


1995; 


Owens, Oosterbroek & Parmar 


2007) in order to investigate the origin of this unusual complex of emission 
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lines. Chandra -HETG observations of 4U 1626-67 made in 2001 revealed a double peaked 

J-II L 

nature of these lines ( Schulz et al . 12001 ) interpreted as Doppler pairs of broadened lines 
of hydrogenic and He like Ne and O. The Doppler velocity corresponding to t he widths of 


these lines is of the same order as Keplarian velocities in the accretion disks. ISchulz et al. 


(2001) have also mentioned the less likely possibility of Doppler shifted lines originating in 
a bipolar disk outflow. In the case of an accretion disk origin, lines of different elements and 
ionisation states will be produced in different annular regions of the accretion disk, at the 
corresponding values of the ionisation parameter. 


Emission lines, especially the iron fluorescence line at 6.4 keV is ubiquituous among the 
X-ray pulsars and has various origins like the acretion disk, strong stellar wind of the high 
mass companion, accretion stream from the companion star, outflowing disk of the Be-star 


companions etc. The line flux is usually found to be variable^ bot 


l with the continuum flux 


and with the abso r ption column den sity fllnoue 


Naik fe Paul 


2003 


Ffirst et al. 


1985 


Makishimal 


1986 


Kotani et ah 


1999 


20111) . However, in the ac cretion powered pulsars, the line 


flux is often found to have no pu 

se pf 

Wilkinson et al. 

2011; 

Suchv et al. 

2012) 


rase dependence (Paul et al. 


2002 


Lei et al.l 


2009 


20121) . An exception is a strong pulse phase dependence 


of the 6.4 keV line seen in Hercules X-l ( Vasco et al. 2013) in which the line flux is found 
to be minimum near the pulse peak and there is a remarkable similarity between the pulse 
phase dependence of the iron line flux and the power-law photon index. 


Pulse phase resolved spectroscopy of 4U 1626-67 carried out in 20 bins show a strong 
flux variation of O VII (0.568 keV) line with pulse phase. We have also shown a comparsion 
of spectrum in two phase ranges (0.0-0.15) and (0.15-0.4) where the flux of this line shows 
maximum difference. However, we would like to add a caveat that if a part of the line of 
sight absorption column is local to the source with a pulse phase dependent ionisation state 
of the material, the same can introduce some unforseen artifacts in the phase resolved mea¬ 
surements of the line parameters. Unfortunately, it is not possible to investigate this with 
a pulse phase resolved spectroscopy of the RGS data. The two lines O VIII (0.653 keV) 
and Ne IX (0.915 keV) show weak variation with phase while the Ne X (1.021 keV) line is 
consistent with being non-varying. 

In 4U 1626-67, we found the O VII line flux to vary by factor of about four (maxi- 
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mum/minimum; Figured] A th panel on the right column), significantly larger compared to 
the relative variation of total flux ( Figure [0 top panel on right column). This indicates 
that the line flux variation is unlikely to be due to varying illumination of flat accretion 
disk. The dimension of the accretion disk in 4U 1626-67, which has an orbital period of 
about 40 minutes is about 0.5 lt-sec which is also smaller than the variability time scale 
observed for the line flux, about 2 seconds. Thus, one possibility of the observed line flux 
variation could be warps in the accretion disk. Milli-Hertz QPOs in the UV and optical 
lightcurves of 4UJ.626-67 have been interpreted as signature of warps in its accretion disk 
(Chakrabarty et ah 20011). In the presence of warps, the visibility of the X-ray illuminated 
part of the accretion disk can have significant variation with pulse phase. If any of the 
emission lines is produced at the same radius as the warp, it may result into a pulse phase 
depen dent line flux. A ve r y low upper limi t of mass function obtained from pulse arrival time 


delay (iL evine et aL 


1988 


Jain et ah 


20081 1 indicates a very small mass of the companion star 


or a large angle of the orbital plane, or a combination of both. While much of the accretion 
disk is likely to be in the orbital plane, the line production region (inner accretion disk) 


probably has warpe d st ructure and 


the mHz QPOs (jChakrabarty et ah 


ras a different inclination, which is also borne out from 


2001 ). 


Though the results presented here are indicative of an accretion disk with a warp around 
a distance at which the O VII emission line is produced, it is also possible to produce a large 
flux variation of one emission line if the accretion disk sees a very different pulse profile 
compared to the observer and if the lines produced have different escape probability. Since, 
only O VII emission line showed strong variation accross the pulse phase, it may be pos¬ 
sible one of the component of triplet O VII emission line may have higher escape probability. 


The X-ray timing characteristics, like the pulse profile, the QPO’s etc of 4U 1626 -67 are 


known to be different i n the spin-down phase compared to spin-up phase flJain. Paul &; Dutta 


2010 


Beri et ah 


2014) which may also be caused by a difference in the inner accretion flow 


from a warped accretion disk in the spin-down phase. If this scenario is true, the Doppler 
characteristics of the O and Ne lines are expected to be different in the current spin-up 
phase of 4U 1626-67. We also note that Her X-l is known to have a warped accretion disk, 
but the pulse phase dependence of the iron line flux in Her X-l is almost identical in the 
different superorbital phases fiVasco et all 20 131 ) which is caused by precession of the warped 
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structure. Thus the variation of the line flux in the two sources 4U 1626-67 and Her X-l 
have different origin and the latter is explained to be due to reprocessing of continuum X-ray 
in the accretion coulmn with hollow cone geometry. 


6 SUMMARY 

We have found the O Vll emission line in the spectrum of the accretion powered pulsar 
4U 1626-67 to vary with pulse phase, the only pulsar after Her X-l to show a pulse dependent 
emission line flux. The flux variation is by a factor of up to 4, much larger than the pulse 
fraction, indicating the presence of warp in the accretion disk at the same radius at which 
this line is produced by reprocessing of the X-rays emitted from the neutron star. This is 
also the first X-ray spectroscopic evidence of an accretion disk warp in an X-ray binary. 
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Figure 1. The left-hand panels show the energy resolved pulse profiles created using EPIC-pn data and binned into 32 
phasebins. The right panel shows the variation of flux (f) of continuum and the low energy emission lines near 1 keV with 
the pulse phase. The units of all continuum fluxes ( ftotali Ibb-, fpL ) are in 10“ 10 ergs cm~ 2 sec -1 , while the line fluxes are 
measured in units of 10 —^photons cm~ 2 sec -1 . All the errors were estimated with 1 <t confidence. 
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Figure 2. Data and folded model of 4U 1626-67 from XMM-Newton EPIC-pn. Upper panel shows the phase averaged 
spectrum. The thick line is the best fit model obtained by fitting the continuum (consisting bbody and power-law). Middle 
panel show the residuals (Counts cm~ 2 s~ 1 keV~ 1 ) obtained, using only continuum models and the lower panel are the residuals 
(Counts cm~ 2 s~ 1 keV~ 1 ) after fitting gaussians at 0.568 keV, 0.653 keV, 0.915 keV and 1.02 keV energies along with continuum 
models. 
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Pulse Phase 


Figure 3. Variation of spectral parameters of continuum accross the pulse phase for 4U 1626-67. The units of Nh (equiv¬ 
alent hydrogen column density) is 10 22 atoms cm ~ 2 . Power-law ( Npl ) and blackbody normalisations ( Nbb ) is in units of 
photons keV -1 cmT 2 s —1 . All the values are in 1 a confidence. 
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Figure 4. Equivalent widths (EW) of Ne and O lines as the function of phase for 4U1626-67. EW is measu 
The errors were estimated with 1 <t confidence. 
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Figure 5. The deviation of the flux of the O VII line with pulse phase with respect to the average flux is shown here. The 
fraction of the number of bins that deviates by more than Ncr is plotted as a function of N. The same for purely statistical 
variation is shown with a line. The y-axis is in log scale. 
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Figure 6. In the left, panel-1 shows data and the best fit models of the spectrum extracted in (0.0-0.15) phase range, the 
middle panel shows the residuals (Counts cm~ 2 sec~ 1 keV~ 1 ) for the best fit and the bottom panel shows the residuals after 
setting line fluxes to zero. Similarly, in the right, the top panel shows data and the best fit models of the spectrum extracted 
in (0.15-0.4) phase range, second panels contains the residuals obtained with best fit while the bottom panel has residuals with 
normalisations set to zero of line fluxes. Note the difference of the flux of the first emission line around 0.569 keV. 
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Figure 7. The variation of flux (f) of O VII emission line accross the pulse phase using different trial periods in 
comparison to the true period=7.67548 seconds. The line flux is measured in units of 10 —^photons cm~ 2 sec -1 . 
All the errors were estimated with 1 a confidence. 




